Reconstruction of interacting dark energy models from parameterizations 
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Models with interacting dark energy can alleviate the cosmic coincidence problem by allowing dark 
matter and dark energy to evolve in a similar fashion. At a fundamental level, these models are 
specified by choosing a functional form for the scalar potential and for the interaction term. However, 
in order to compare to observational data it is usually more convenient to use parameterizations 
of the dark energy equation of state and the evolution of the dark matter energy density. Once 
the relevant parameters are fitted it is important to obtain the shape of the fundamental functions. 
In this paper I show how to reconstruct the scalar potential and the scalar interaction with dark 
matter from general parameterizations. I give a few examples and show that it is possible for the 
effective equation of state for the scalar field to cross the phantom barrier when interactions are 
allowed. I analyze the uncertainties in the reconstructed potential arising from foreseen errors in 
the estimation of fit parameters and point out that a Yukawa-like linear interaction results from a 
simple parameterization of the coupling. 

PACS numbers: 98.80.Cq 



INTRODUCTION 

Recent data from type la supernovae (SNIa) [l], Q, 
cosmic microwave background (CMB) [3j and large scale 
structure (LSS) Q all point to the fact that the universe 
has recently entered into a stage of accelerated expansion. 
This is an unexpected and revolutionary discovery that 
calls for new physics, since gravity is always an attractive 
force. 

The simplest possibility to explain the acceleration of 
the universe is to postulate the existence of a cosmo- 
logical constant of the right magnitude. While this as- 
sumption is still compatible with all available data, it is 
unsatisfactory on the grounds of the huge amount of fine 
tuning required. Hence, exploratory models with scalar 
fields possessing a varying non-zero vacuum energy den- 
sity, usually called quintessence models have been 
studied as an alternative to the cosmological constant 
solution. One of the most important tasks ahead of ob- 
servational cosmology is to devise methods and gather 
data to successfully and unequivocally distinguish be- 
tween these two possibilities. Finding evidence for an 
evolving vacuum energy would be one of the greatest dis- 
coveries of the century. 

What is the appropriate scalar potential that could 
reproduce the observed data? There has been a large 
number of papers dealing with the possibility of recon- 
structing the potential directly from the data @, 0, H, 
0, 0, Ell 0, El, E3 and it is fair to say that any state- 
ment about the form of the potential depends upon some 
level of parameterization. The procedure has also been 
extended to reconstruct the potential and the form of the 
scalar-gravity coupling in the Jordan frame of a scalar- 
tensor gravity theory from data on luminosity distance 
and linear density perturbation (l5j . More recently Guo, 
Ohta and Zhang [l6( studied the reconstruction of the 



scalar potential directly from parameterizations of the 
equation of state. 

An intriguing possibility is that the scalar field is not 
totally decoupled from our world [13] • Although the cou- 
pling to baryonic matter is severely constrained, it is still 
possible to allow for a coupling to non-baryonic dark mat- 
ter. From a lagrangian point of view, these couplings 
could be of the form W (fymQ^ip or W(<fi)mQtp 2 for a 
fcrmionic or bosonic dark matter represented by ip and 
<p respectively, where the function W of the quintessence 
field 4> can in principle be arbitrary. In this scenario, the 
mass of the dark matter particles evolves according to 
some function of the dark energy field <f>, leading to an 
effective equation-of-state for the dark matter. 

I will consider exploratory models for dark energy with 
a canonical scalar field coupled to dark matter. The la- 
grangian for this class of models is specified by the choice 
of two functions of the scalar field: the scalar field poten- 
tial V((f>) and W(<p), the function that characterizes the 
coupling to dark matter. Specific forms for the scalar 
potential and for the interaction, such as power-law or 
exponential functions, have been extensively studied in 



the literature HI, 



However, instead of postulating a concrete model by 
choosing definite parametric forms for V(<fi) and W (</>), it 
is often more convenient and completely equivalent to in- 
troduce a time-dependent parameterization for the dark 
energy equation of state wde{o) and for a coupling func- 
tion 5(a), where a(t) is the scale factor of the universe. 
In fact, this approach is widely used in the uncoupled 
case in order to test for the time- variation of dark energy 
since it is much easier to compare to SNIa observations 



2l| . For the coupled case, an analysis of the impact of 



interactions on SNIa observations was performed using 
a sp ecific coupling either in the Einstein [221 ] or Jordan 
(231 ] frame of a scalar-tensor theory. More recently, pa- 



2 



rameterizations of wde(o) and 6(a) were directly used 
to study the effects of DE-DM interactions on fits from 
SNla data Q. 

Given a parameterization for wde^) and 6(a) with 
parameters fitted by observations, it is important to in- 
vestigate the shape of the potential and interaction func- 
tions that results in those parameterizations. They define 
the fundamental model behind the parameterizations. In 
this article 1 will show how to reconstruct the scalar po- 
tential and the interaction from the general parameter- 
ized forms of the dark energy equation of state wde(o) 
and a coupling function 6(a) and provide a few concrete 
examples. 



RECONSTRUCTION 

The reconstruction program proposed here generalizes 
the one first developed by Ellis and Madsen [25| . I con- 
sider a spatially flat universe composed of three perfect 
fluids, namely dark energy, non-baryonic dark matter 
and baryons. The dark matter and baryons are non- 
rclativistic pressureless fluids and Einstein's equations 
result in: 



H z = 



8ttG 



(P4> + RDM + Pb) 
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H + H 2 = — (p^ + p DM + Pb + 3pct>) 



(1) 



where H = a/a and the dark energy is described by a 
scalar field (f> with energy density and pressure given by: 

P^ = \tf + V(&, p^^-V^) (2) 

with a equation of state w^, defined by — p^/p^. 

Conservation of the stress-energy tensor requires that 
the total energy density and pressure obey: 



pT + 3H(p T + Pt) = 0. 



(3) 



Introducing the coupling function 6(a) between dark en- 
ergy and dark matter as 



6(a) 



d In nil/, (a) 
din a 



(4) 



results in the following equation for the evolution of the 
DM energy density p^M [261 ]: 

p DM + iHp DM - 6(a)Hp DM = 0. (5) 

Also, conservation of baryon number requires: 

p b + 'AH Pb = (6) 

and eq. ([3]) then implies that the dark energy density 
should obey 



Notice that the parameterization eq. ([4} implies 

W(cj>(a)) = e-Sl &{a ' )dXna ' (8) 

normalized such that W(cf>(a = 1)) = 1. 

Combining eqs. (|2l7l8p one obtains a modified Klein- 
Gordon equation for the scalar field: 

^+(f + *^)=». « 

in agreement with Das, Corasaniti and Khoury jl9j |. 

One can now proceed to reconstruct the potential and 
the interaction for a given parameterization of the equa- 
tion of state w(a) and the interaction 6(a). The first 
step is to find the time variation of dark matter energy 
density, which is easily obtained by solving eq. (|S"|): 

( \ (°) -3 - f 1 S(a')d\aa' /in\ 

Pdu(o) = p y D ' M a J e J» , (10) 

where p^ M is the non-baryonic DM energy density today. 
It is more useful to work with the variable u — In a, and 
one can write 



Pdm 



(u)=p { ° ) M e~^e-f>^ du '. (11) 



The second step is to substitute pdm(u) into eq. ([7]), 
which in terms of u reads: 

p'^u) + 3(1 + w^(u)) P<l> (u) + 5(u) PDM (u) = 0, (12) 

where ' = d/du, and find a solution p^(u) with initial 
condition p$(u = 0) = p^\ with /A ' being the dark 
energy density today. 

In the third step one constructs the Hubble parameter: 



H 2 (u) 



-3„ e - + n ,/( u ), 



(13) 

where fix = pjP / p^ , the critical density today is p^ = 
3Hq / (8irG) and Hq is the Hubble constant. The function 
f(u) that determines the evolution of the dark energy 
density is in general obtained numerically. 

Having obtained the Hubble parameter, the fourth step 
consists in solving the evolution equation for the scalar 
field obtained from eqs. fH 




1 / d\nH(u) 3 , . _ . 

' y ^ + -(n DM (u) + n b (u)) 

(14) 

where </> = tfi/Mpi is the scalar field in units of the Planck 
mass (Mpi = l/-\/G) and 



P0 + 3ff(p0 + P4) + 6(a)Hp DM = 0. 



(7) 



&DM,b(u) 



PDM,b(u) 



P<t>(u) + Pdm(u) + Pb(u) ' 



(15) 



3 




FIG. 1: Reconstructed potential for — —0.99 and 8 = 0. 




In the fifth step one numerically inverts the solution <p(u) 
in order to determine u{(j)) to finally obtain 



V{4>) 



and 



o (0) 



1 H(u) dH/H 
3 Hq du 



H 2 {u) 



-il b e - -ilDAie e 



w(u(4>)) 



i>0 . 



(16) 

J S(u')du' 

(17) 



This completes the reconstruction procedure. I will now 
work out some examples of this procedure. I adopt r2^ = 
0.7, £Idm = 0.25 and Oj, = 0.05 in the following. In 
all examples I integrate the field equation starting from 
Ui = —1.8, corresponding to Zi = 5.05 and I arbitrarily 
set <f>(ui) = —1. 



EXAMPLES 
Constant w^, and 5 

I start by considering the simple example of a constant 
equation of state and constant coupling d. In order to 
gain some intuition, I first show in fig.[T]the reconstructed 
potential for wq = —0.99 and 6 = 0. As expected, the 
scalar field evolves slowly in a very flat potential. Notice 
that the potential energy is slightly below 0.7p c today 
due to the small kinetic energy of the scalar field. 

I now turn on the interaction. In this case one has 301: 



/ \ (0) -3+5 

Pdm {pi) = p DM a ^ , 



(18) 



and the solution to eq. (|7j) is 

p,{a)=p^a-^+^ + 



S (°) („-3(l+v>4,) _ „-3+<5 



S + 3ui0 



Pdm a 



FIG. 2: Reconstructed potential for W4, = —0.9 and 5 
(solid), S = 0.05 (dashed) and 5 = 0.1 (dot-dashed). 



The first term of the solution is the usual evolution of 
DE without the coupling to DM. From this solution it 
is easy to see that one must require a positive value of 
the coupling i5 > in order to have a consistent positive 
value of for earlier epochs of the universe. This feature 
remains in the case of varying and in the rest of the 
paper I will assume that S is positive. 

In fig. [2] I show the effects of the coupling in the re- 
constructed potential. I use wo = —0.9 and 5 = 0, 0.05 
and 0.1. The scalar potential becomes steeper for larger 
values of the coupling due to the different dynamics in- 
troduced by the DE-DM coupling. 

One might have expected that the introduction of the 
coupling could allow the possibility of having a phantom 
equation of state, < — 1. However, it is easy to show 
that one can write eq. (fT4"|) for both the uncoupled and 
coupled cases as: 



du 



■ w<f,(u)) n<j,(u) 



(20) 



in the general case of a time-varying equation of state 
and hence, in order to have a real scalar field one must 
consider only w^{u) > — 1. 

However, the fitter who is unaware of the interaction 
would instead find an effective equation of state w e f / (u) 
defined implicitly by: 

p (0) e - 3( l +WefAu))u + p (0) Me - 3u = ^ + p (0) Me (S + S)u_ 

(21) 

In fig. [3] I show that the effective equation of state can 
in fact cross the phantom barrier. This possibility was 
also recently pointed out in the context of scalar-tensor 
theories of gravity [27j ]. 

One can also easily reconstruct the interaction W(<p) 
in this simple case: 



(19) 



w(4>(u)) 



„5u 



(22) 
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FIG. 3: Effective equation of state for w = —0.9 and 8 = 0.1. 
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FIG. 4: Reconstructed interaction for w = —0.9 and 8 — 0.05 
(dashed line), 8 — 0.1 (dot-dashed line) and w — —0.99 and 
8 = 0.1 (solid line). 

which of course guarantees that W(<p(u — 0)) = 1 as 
required. In fig. [5] I plot the reconstructed interaction 
term for w = —0.9 and S = 0.05, 0.1 and also for w = 
—0.99 and 5 = 0.1. Notice that the interaction increases 
with time, corresponding to a mass that decreases with 
increasing redshift. One can see that the interaction term 
becomes steeper with increasing 5 and decreasing w. 



FIG. 5: Reconstructed potential for u>o = —0.91, w\ = 0.4 
and 5 = (solid line). The allowed region from the 1 — a 
uncertainties in the estimation of the equation of state pa- 
rameters is the region between the dashed lines. 
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FIG. 6: Reconstructed potential for u>o = —0.91, w\ = 0.4 
and 8 — 0.1 (solid line). The allowed region from the 1 — a 
uncertainties in the estimation of the equation of state pa- 
rameters is the region between the dashed lines. 



Variable and constant 8 

I now analyze the potential reconstructed from the 
commonly used 2-parameter description of the equation 
of state [28j |: 

w<t,{a) = wo + 101,(1 — a). (23) 

A recent fit to SNIa, CMB and LSS data obtained (with- 
out DE perturbation) [iij ]: 

wo = -1.098i°;°8o, tyi = 0.416±8:?H (24) 

Hence I will assume that it is possible to fit the 
2-parameter equation of state with a 1-a precision of 



6wo = 0.08 and Swi = 0.20. For illustration purposes, 
I will take the central values wq = —0.91 and w\ = 0.4 
and study the effects of dark energy interaction in the 
reconstruction of the potential. In fig. [5] I show the al- 
lowed region in the potential form arising from the 1-a 
uncertainties in the estimation of the equation of state 
parameters. This region therefore gives an idea of the 
uncertainty one can expect in the reconstructed poten- 
tial given the foreseen errors in the fit parameters. 

In fig. [6] I show the effect of interaction with 6 = 0.1 
on the allowed region in the potential form arising from 
the 1 — a uncertainties in the estimation of the equation 
of state parameters. 
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-0.98 -0.96 -0.94 -0.92 

FIG. 7: Reconstructed potential for constant w = —0.9, and 
So = 0.1 (solid line) and So = 0.2 (dashed line). 

Constant w^, and variable S 

Finally, I analyze the case of a constant equation of 
state and an interaction parameterized by the function 

8(a) =d -^, (25) 
1 + or 

which is well behaved in the past as well as in the future 
and 5 (a = 1) = 5q. In this case one finds: 

PDM (a) = p ^ a -3 e 25o(axct a n(a)-./4) (2g) 

The reconstructed potential is shown in fig. [7] for = 
—0.9 and 6q = 0.1,0.2. The interaction function is given 
by 

W(u(4>)) = e 2«o(arctan(e"W))- W /4) ^ 

and it is plotted in fig. [H Notice that in this example 
the interaction function has an approximate Yukawa-like 
linear form, even more so for the 5q = 0.2 case. 

CONCLUSIONS 

It is natural to study models where the dark energy 
field is not totally separated from the rest of the world. 
In principle, it can interact with dark matter and this 
possibility has interesting consequences in the evolution 
of the universe, such as a mass- varying cold dark matter 
particle with a non-zero effective equation of state. 

In order to compare these models to observations it is 
more convenient to use parameterizations of the dark en- 
ergy equation of state and the coupling to dark matter. 
However, once the parameters are estimated it is impor- 
tant to find the fundamental lagrangian of the theory, 
that is, to determine the functional form of the scalar 
potential and its interaction with dark matter. 




4> 



FIG. 8: Reconstructed interaction function for constant w = 
—0.9, and So = 0.1 (solid line) and <5o = 0.2 (dashed line). 



I showed in this paper how to reconstruct the scalar 
potential and the scalar interaction with dark matter 
from general parameterizations. For illustration purposes 
some examples are worked out. Uncertainties in the re- 
construction due to uncertainties in the estimation of pa- 
rameters are analyzed. It is pointed out that the phan- 
tom barrier can be crossed if the fit does not take into 
account interactions and that a Yukawa-like linear in- 
teraction results from a simple parameterization of the 
coupling. 

More precise data from CMB, SNIa and LSS from dif- 
ferent collaborations is expected to arrive in the near 
future. These data can be used to estimate the parame- 
ters of simple parameterizations and the procedure shown 
here can at least point towards the underlying fundamen- 
tal model describing our universe. 
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